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Abstract 

The paper considers pseudo-differential boundary value control systems. The underlying 
operators form an algebra T> with the help of which we are able to formulate typical boundary 
value control problems. The symbolic calculus gives tools to form e.g. compositions, formal 
adjoints, generalized right or left inverses and compatibility conditions. By a parametrizability 
we mean that for a given control system Au = one finds an operator S such that Au = 
if and only if u — Sf. The computation rules of T> (or its appropriate subalgebra T>') 
guarantee that in many applications S can be refinery analyzed or even explicitly calculated. 
We outline some methods of homological algebra for the study of parametrization S. Especially 
the projectivity of a certain factor module (defined by the system equations) implies the 
parametrizability. We give some examples to illustrate our computational methods. 



AMS-classification. 93B25, 93C20, 93A10, 35S99 

1 Introduction 

The general theory of control systems corresponding to boundary value problems for linear and 
nonlinear partial differential equations (PDEs) has been developed by two apparently different 
approaches: One has applied functional analytic and algebraic methods. In addition, differential 
geometry increasingly offers effective tools and gives geometric intuition in this field. Functional 
analytic methods are based on the analysis of underlying operators defined abstractly in appropriate 
(Banach) spaces, e.g. [TJ [4j [3] . Definitions of basic concepts such as controllability, observability, 
stability are typically formulated using trajectories which have direct intuitions to the real world 
problems in technology and in science. 

In frequence space problems (e.g. transfer function analysis) algebraic approach is more con- 
ventional [231 [21] but in state space problems algebraic study of control systems is quite recent 
method. The algebraic methods analyze algebraic structures such as structures of certain modules, 
exact sequences, compatibility conditions and one-sided inverses. In these approaches one mixes 
the input, state and output variables (which are called control variables) and studies only the 
structural properties of this mixed system (e.g. [2H [201 [3] ) ■ In the certain linear cases one is 
able e.g. to apply the related (differential) algebraic modules to determine whether the system is 
controllable or not. The algebraic methods have often less intuitive connections to real problems. 

Here we consider parametrizability of certain linear control systems which means that for a 
given control system Au = one finds an operator S such that Au = if and only if u = Sf. The 
operators A and S are built on the operators of V where V is the chosen algebra of operators. 
Parametrizability is a kind of structural (internal) property of the control system. The concept 

'Department of Mathematics and Statistics, University of Kuopio, P.O.Box 1627, FIN-70211 Kuopio, Finland. 
Email: Jouko.Tervo@uku.fi 



1 



is related to the differential flatness of (ordinary differential) control systems ( [21 E] [T2] ) . The 
"parameter /" is in the role of flat output. / is endogeneous if it can be conversely expressed with 
the help of u by / = Qu where Q is built on the operators of V. In this case the system can be 
considered to be " flat" . In some cases parametrizability is equivalent to the torsion freeness of 
certain structural modules (|20jV 

By now the structural study of state space problems, especially from the algebraic point of 
view, is mainly developed without boundary values. Some attempts nevertheless can be found in 
literature [6l HH [15] . It is clear that the structural properties related to the boundary value control 
problems depend also on the boundary conditions, not only the PDE system. In the following we 
suggest how to take into account the boundary conditions in parametrization. We use very general 
formulation of boundary value systems. The operators consists of appropriate pseudo-differential 
and boundary value operators. This calculus enables, for example, the consideration of the trace 
operator. In addition, the calculus gives the greater freedom in the manipulation of e.g. the 
compositions, adjoints, compatbility conditions and left/right inverses. We express our formalism 
tailored for the boundary value problems related to certain partial differential equations. Besides 
boundary value control systems for PDEs, our approach enables also e.g. the consideration of delay 
systems. 

Firstly we choose an appropriate algebra T> of pseudo-differential and boundary value operators. 
This algebra originated in [14] and afterwards it has been enlarged ([3 [9]). The algebra gives very 
natural frame to calculate with PDEs and related boundary value operators. We generalize the 
definition of intrinsic parametrization (the word "intrinsic" refers to the situation where state, 
input and output variables are mixed) for the control systems which are corresponding to this 
algebra and we give some preliminary tools to study and find the parametrization. Compared 
with the earlier differential algebraic methods we loose certain structures such as differential fields 
and the derivation rules typical in the context of partial differential operators. In our case the 
modules are not generally over entire rings (integral domains). This follows e.g. because of the 
ring C°°(G x A) is not entire: We have V\V2 = for any nonzero functions whose supports are 
disjoint. Furthermore T> as the ring of matrices is not entire. 

In section 4 we formulate some preliminary ideas to study the parametrizability by certain 
homological algebraic methods. We define (more general) D'-parametrizability concept, where 
£>' is a subalgebra with unity of T>. We don't try to make clear here relations between (various 
variants of) controllability and parametrizability. The definition of structural controllability for 
the boundary value control systems (under consideration) can be potentially founded on the use 
of torsion freeness concept which can be generalized also for nonintegral domains. Because the 
underlying module is not generally over an entire ring the application of torsion concept is slightly 
more complicated than in the conventional case. 

One aim of this paper is to show that the applied algebra gives tools to construct the operators 
needed in parametrization. Parametrizabilty gives potential methods for practical controller design 
since from u = Sf it follows that we can express, say the input variables c = (c\, ...,c p ), the state 
variables v = (vx, v q ) and the output variables y = (yi, ...,y r ) with the help of " (free) variables 
/" as 

c=Sr/, v = S 2 f, y = S 3 f. (1) 

Here the operators Sj are constitued of the operators lying in the algebra V and so we can calculate 
c, v, y using computational rules of T>. For example, in the output tracking problem one seeks the 
"parameter" / = /* such that (at least optimally) y* = S3/* where y* is a given reference output. 
The required input is simply c — c* = Si/*. This kind of "open controls" are important in some 
problems of modern technology (in molecular physics, for example) . The other potential field of 
applications is the optimal control problems. Roughly speaking this option can be described as 
follows: If we have an optimal control problem min u F{u) under the contraint Ait = 0, we can 
transform it to unconstrained problem min/ F(Sf). 
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1.1 Basic notations 

Let G be an open bounded set in R n and let A be an interval in R. We assume that the closure 
G is a smooth differentiable manifold with boundary. Furthermore we suppose that the boundary 
dG is orientable and that the unit normal vector v = v(x), x G dG is pointing "outwardly" on 

dG - _ 

The spaces C°°(G), G°°(<9G), G°°(G x_A) and C°°(dG x A) are correspondingly the collections 
of smooth functions G — > R, dG — > R, G x A — » R and dG xA^R. C£°(G) is the space of 
test functions and S(R n ) denotes the Schwartz space of rapidly decreasing functions R" — > R. 
Furthermore, we define 

G°°(G x R+) = {(f) G C*°°(G x R+) | for which there 

exists a ip £ G°°(G x R) such that = V'Igxr } 

and 

C£°(G x R + ) = {0 G C°°(G x R + ) | has a compact support}. 

D'(G) is the space of distributions that is, D'(G) is the dual of G£°(G). 
Denote 

In addition, we denote 

ft = ft 1 ' 1 = G°°(G x A) © G°°(o>G x A). 

Let 

j^N,m _ {( VfW ^ g 7?. Ar ' m | u and io have compact supports 
in G x A and in dG x A, respectively}. 

For (v,w) G ^ JV ' m , (<f),ip) G ft^'™ we can define an inner product 

N m 

({v,w),(4>,il>)) = ^2(vj,<f>j) L2(GxA) ( w k^k) L2(eGxA) ■ (2) 

j=l k=l 

Let u (or Tu) denote the standard Fourier transform 

&(0 = = / u^Je-^^da! (3) 

for a function u in the Lebesgue space Li(R"). The inverse Fourier transform is for appropriate 
functions 

The ring of m x n-matrices is denoted by Af(m x n). The partial derivative is denoted more 
shortly by d" . 

2 Control system 

2.1 Pseudo-differential and boundary value operators 

We formulate the (generalized) boundary control systems with the help of pseudo-differential and 
boundary value operators (e.g. [BE]). Let S™(G x A x R n+1 ) be a space of functions, so called 
anisotropic symbols, a G G°°(G x A) such that for any compact set K C G x A, a, j3 G Nq, j, I G 
No there exists a constant C^pj^jc > such that 

\d^didld l v a{x,t,i,ri)\ < C at0 , jtltK p(Z, V ) m -W- vl for all (£,77) G R" +1 , (a:,t) G X (4) 
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where 

p(tv) = (i + \t;\ 2 + \v\ 2/ n 1/2 - 

As is standard, the space of symbols S™(G x A x R™ +1 ) can be equipped with the topology 
defined by the appropriate semi-norms. In the case where we have no boundary considerations, 
the (classical) pseudo-differential operator A : C^°(G x A) — > G°°(G x A) is defined by 

Af>{x,t) = (2n)-w / / a{xM,v)kt,v)*£<*n- (5) 

JR" JR 

When the boundary dG is included in the considerations the definition of pseudo-differential op- 
erator contains some modifications. Especially, the so called transmission condition is needed. In 
addition, we need several classes of other operators (boundary operators). 

It is sufficient to explain the operator classes only in the local case that is, in the case when 
the closure G is replaced with a set G' x R + where G' is an open subset of R" -1 and where the 
boundary dG is replaced with a set G'. The technically tedious reduction to that case is based on 
the partition of unity and it is standard. 

Generally speaking, the symbols of "boundary operators" satisfy conditions analogous, in a 
sense, to the transmission condition. This provides the inclusion 

A(C™(G' x R + x A) © C C °°(G' x A)) c C°°(G' x R+ x A) © G°°(G' x A) 

for any boundary value operator A explained in section f2. 2. 11 In addition, the kernels of operators 
involved in a boundary value problem must have proper supports. Then 

A(C°°{G' x R + x A) © C°°(G' x A)) c C°°{G' x R+ x A) © G°°(G' x A) 

and one can compose the operators of boundary problems. In the following we explain the operators 
in detail. 

A pseudodifferential operator (ipdo) 

r + A: G C °°(G' x R+ x A) -> D'{G' x R+ x A) 

is defined by 

r+Av{x,t) = (27r)-(" +1 V+ f [ e l{{x ' t) ^^) a (x,t,Z,T 1 )(e~+v){£,r 1 )dt;dT 1 . 

JR" Jr 

The symbol a for any compact K C G' x R + x A satisfies 

\d^dtd^d l n a(x,t,C,v)\ < C a , jW p(tv) m - m ~ vl , (xM,v) £Kx R"+\ (6) 

where a, (3 G Nq , j, I G N . The space of symbols satisfying §6§ is denoted by S™(G' x R + x A x 
R™ +1 ). The number m is called the order of the symbol a and, simultaneously, the order of ^do 
r + A. Above r + refers to the restriction operator r + f = ,/|g'xr + xA and e + refers to the extension 
by zero from G' x R + x A on G' x R x A . 

The symbol a satisfies the transmission condition, if the expansion 

d% n a{x,t,£,r])\ Xn=0 = '^2a 1P {x',t,t',r))tP+ ^ a jk (x' , t, rj) — ^ — (7) 
P =o k=-oo (PU'iW+^n) 

holds for all 7 G No while a lv G S™~ P (G' x A x R"), a 7 fc is rapidly decreasing sequence in 
S'™ +1 (G' x A x R") (that is, for any semi-norm p on 5™ +1 (G' x A x R") and any N G N there 
exists a constant c N , P such that p{a lk ) < c N . p (l + k)~ N ) and p{g ,T}) = (1+ |£'| 2 + \n\ 2 / u ) 1 / 2 . The 
space of symbols of order m satisfying the transmission condition will be denoted by U" 1 . 
Assertion 1. If a tjjdo r + A satisfies the transmission condition then r + Av G G°°(G' x R + x A) 
for any v G G C °°(G / x R+ x A). 
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The consideration of Assertion 1 can be reduced to the one-dimensional case in the same way 
as in [14], Theorem 2.9. Partial differential operators are typical examples of pseudo-differentials 
operators. 

A potential operator. A function k S C°°(G" x R + x A x R" +1 ) is called a potential symbol 
of order m if 



k(x, t, (;', £„, v) = £ a k (x, t, V ) (p(£, r,) - iQ" (p(£', »?) + , (8) 



where dfe is a rapidly decreasing sequence in the space S 7 J l+1 (G' x R + x A x R") consisting of 
functions subjected to © with 77) replaced by rf). The space of potential symbols of order 
m is denoted by K™. 

A potential operator is defined by 

Kw(x, t) = (2^" +1 > / / e^ x ^^k(x, t, H\ £„, »7)t&te', »?) ^ 

for w e C C °°(G' x A). It is obvious that if(C c °°(G" x A)) c C°°(G' x R+ x A). The (classical) 
solution operators of certain partial differential boundary value problems contain integral terms 

Kg(x,t)= [ G(x,t,y',T)g(y',r)dadT (9) 



which are examples of potential operators. The other examples are the adjoints of the below 
defined trace operators which are often potential operators. 

A trace operator. A function t £ C°°(G' x A x R™ +1 ) is called a trace symbol of order m 
and class d if 

t(x', t, £„, r,) = «p(^, t, rj)% + £ a ^'> *' £'> r?) (p( ^ ,7?)+ ^ +1 (10) 

where a p 6 S™~ P (G' x A x R") and afc is a rapidly decreasing sequence in the space 5™ +1 (G" x 
A x R n ). The space of trace symbols of order m and class d is denoted by T™ ,d . 
A trace operator is defined by 

Tv{x' ,t) = {2«)- n f f e *<(-'.*).«'.")> [ t(x\t,e,Z n ,r,)(e^v)(Z, V )dt n dt'dr, 

JR"- 1 JR J 

for t> e [G' x R + x A) . Here f + fd£ n is the curve integral f f(z)dz where 7 is " a path rounding 

the upper complex half plane Imz > in the counterclockwise direction" that is, f + fd^ n = 
liniR^oo § Tr f(z)dz where Tr is the boundary of the half ball -8(0, R) fl {z € C| Imz > 0}. The 

operator T can be written in the form Tv = X)p=o ^p(px n v { x ' > 0) 0) + Tqv where S p are ipdos with 
symbols in S™- p {G' x A x R") and T is of class 0. It is obvious that T(C C °°(G' x R + x A)) C 
C°°(G' x A). Restrictions of partial differential operators on the boundary (which usually appear 
in boundary conditions) are typical trace operators. 

A Green operator. A function b 6 C°°(G' x R + x A x R™+ 2 ) is called a Green symbol of 
order m and class d if 



b(x,t,e,tinXn,V) (11) 



l+l ' 



where k p £ K™ p and dji is a rapidly decreasing double sequence in the space S™ +2 {G' x A x R n ). 
The space of Green symbols of order m and class d is denoted by B™ ,<z . 
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A Green operator is denned by 

Bv(x,t) = (27T)-"- 1 / / e *«^).(^)> / b(x,t,e,Z n ,( n , V )(^u)(e,( n ,r))d( n cl4dr, 



for v e C c °°(G'xR + x A). The operator £ can be written in the form Bv = K p{ d % n v ( x '> °> t )) + 

Bqv where K p are potential operators with symbols in K™ _p and Bq is of class 0. It is obvious 
that B{C^(G' x R + x A)) C C°°(G" x R + x A). Typical sources of singular Green operators 
are the compositions of potential and trace operators. They are also born in the compositions 
of truncated pseudo-differential operators. Finally they appear e.g. in many solution operators 
related to the usual boundary value problems. 

A pseudodifferential operator on the boundary 

Q: C C °°(G" x A) -> C°°(G" x A) 

is defined by 

Qw(x',t) = (2tt)-" / / e^'^^'^qix'^^^^wiCvWdri. 

JR n ~ 1 Jr 

with the symbol q £ S™(G' x A x R n ). Partial differential operators on the boundary (manifold) 
form an example of these operators. 

In practice the operator r + A is often a partial differential operator. The operator T is a 
natural generalization of the usual partial differential trace operator appeared in classical theory 
of boundary value problems. The operators K, B and Q are needed e.g. in the consideration of 
the compositions, inverses and adjoints. The operators T, K, B, Q are called boundary operators. 
Together with these operators one can formulate very rich variety of boundary value problems. This 
more general formulation also gives more symmetry in the calculus of adjoints and compositions. 

Example 1. Let L(D) — ad x + bd x + c be a PDO with constant coefficients and let d%, d 2 be 
constants (in this example we have no time variable). Furthermore, let G =]0, 1[c R. Then under 
relevant assumptions the solution V\ for the system 

L(D) Vl = q(x)v 2 (12) 
d xVl (0) + d lVl (0) = 0, x «i(l) + d 2Vl (l) = (13) 

is formally given by (U5]j ) 

v = P (qv 2 ) + R (qv 2 ) (14) 
where Pq is a pseudo- differential operator 

P^x) = (27T)- 1 j e*<*.0 jTxfa)<% (15) 

and Rq is a singular Green operator (note that in one- dimensional case the integration reduces to 
summation) 

R ° m = -^^4+fla(o.o)W(i)-5a(o,iW(o)] 

~l^[9i(0,0)W(l) - (~ +5i(0,l)W(0)]. (16) 
Here L(£) — — a£ 2 + M£ — c is the symbol of L(D) and 



9k(x, A) = (2tt) *pv 



X-L(O a \ 



, A = 0,1, k = l,2, 



IR 

W(0) = 51 (0, 0)52(1, 0) - (1 + 5i (1, 0))(1 + 52 (0, 0)), 



where 



9i(€) = ( ad 2 - b- ai£)e g 2 (£) = b - ad x + ai£. 
Here pv denotes that the integral is taken in the sense of principal value. 
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2.2 Computation rules 

We survey some basic computational properties of the above defined operators. 



2.2.1 Operator algebra 

fr + A + B K\ 

A boundary value operator A — \ ^ q J is called proper if the kernels of all operators 

r + A, B, K, T, Q have proper supports. The space of proper boundary value operators is denoted 
byV 

Assertion 2. If A and B are in V, then AB := A o B G V. 

To verify this assertion one can modify the arguments used in [7] for the case of parameter 
dependent boundary value problems. 

(r + A + B K\ 
^ q J is of class V" l ' d , 



m = 



mi, m 2 m 3 \ _ ( d 2 



TO4 ?7l5 I ' \d,4 



when 

the symbol of r + A is in the class U™ 1 , 
the symbol of B is in the class B™ 2,d2 , 
the symbol of K is in the class K™ 3 , 
the symbol of T is in the class T™ 4 ' d4 , 
the symbol of Q is in the class S"™ 5 . 

The assertion 2 implies that proper boundary value operators of arbitrary order m and class d 
form an algebra T> with respect to standard addition and composition of operators. In addition, 
from the above definitions and assertions it follows that the spaces V and C°°(G' x R + x A) © 
C°°(G' x A) are V-modules. The order and the class of AB can be calculated with the help of 
orders and classes of A and B. 

In the global setting we have operators 

r+A : C°°(G x A) -> C°°(G x A) 

K : C°°{dG x A) -> C°°(G x A) 
T : C°°{G x A) -> C°°(<9G x A) 
B : C°°(G x A) -> C°°(G x A) 
Q : G°°(9G x A) — > G°°(<9G x A) 

and 

.4 : ft -> ?e 

where 72. = G°°(G x A) G°°(aG x A). In addition TZ is a X>-module. 
2.2.2 Formal adjoint 

As mentioned above in 7?. = G°° (G x A) © G°° (3G x A) we use the inner product (e.g. for functions 
Ux eTZ, U 2 e TZ Q ) 

(Ui,U 2 ) = {vi,v 2 ) L2(GxA) + (w 1 ,w 2 ) L2(dGxA) 

where 

U k = (v k ,Wk), k= 1,2. 
For any boundary value operator A G T> there exists an operator 

B : G C °°(G x A) G C °°(<9G x A) -> D'(G x A) © D'{dG x A) 
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such that (AU, V) = (U, BV) where (•, •) is the standard duality form on 

(C* C °°(G x A) © C™(dG x A)) x (D'(G x A) © D'(dG x A)). 

If the operator B maps the space 1Z into itself we will say that formal adjoint to A exists and will 
denote the operator B by A* . When the operator A* : 7Z — > 7?. exists we have 

(y,^c/) = (^i*v,f/) (17) 

for V,U G 1Z . 

Assertion 3. The adjoint operator A* exists and is in the algebra T> if and only if the operator 
fr + A + B K\ 

A = I „ n I is of class (i.e., the class of operators B and T equals 0) and in the 



T Q, 

transmission condition |?]] for the symbol of r + A all a lp , except a 7 o, are equal to 0. In the case 
where A* exists then there exists A** and A** = A. 

We omit here the proof of Assertion 3 but we remark that analogous results and considerations 
can be found in the monograph [7]. 

We find that for an operator A G D the formal adjoint A*, generally speaking, does not exist. 
However, one can consider subset T>\ C T> that consists of boundary value operators satisfying 
conditions of the Assertion 3. Then T>\ is an algebra and for any A G T>\ there exists the formal 

adjoint A* in T) x . If wc fix m Q = and d ° = (o) ' thcn the subal S ebra V ™°' d ° c P i 

possesses the same properties as T>\, 

Example 2. Let G =]0, l[c R and 

Tv = d(d)v\ dG 
where d(d) is a pseudo-differential operator 



R 



d(d)v(x) = (27T)- 1 / e^ x d(0e+v(0d^ 
such that its symbol d(£) satisfies 



\d(0\<C-^^, k>1. (18) 



Then we have for v G C°°(G), g G C°°(dG) 



( t ^9)l 2 (og) = (Tv){0)g{0) + (Tv)(l)g(l) 



(2^)- 1 .g(0) / d(0e+t^)d£+ (270-^(1) / e^d&e+viO^ 
Jr Jr 

p 1 ^ / ^ 7777TT / - '-'.•/(< wi.< )n i,i >■ 4- / i'TTTTT / ,^--> + i>, 



(27T)" 1 (JU(0) J^e-^d((;)dt)v(x)dx + J[g(l) J^- x+1 >d(^ ) v(x)dx 



where 



T*g(x) = (27T)" 1 (^g(0) J^ x d(0dC + g(l) J R e^- x+1 U(0d^ 



(19) 



We see that T* is a potential-type operator, but it does not generally possess the transmission 
property. Consideration of operators like H19\) can be found in \25T/ . 
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2.2.3 Order reducing operators 

The so called order reduction can be used to help the forming of adjoints in parametrization 
processes. 

Let A be a boundary value operator of order m = [ TOl ' 771-2 777-3 | and class d = | °! 2 
If T = ^ r+ Q Tl y J while ordr+Ti = -Ni (that is, the order of r+Ti is -JVi) and ordT 2 = 
-N 2 , Nu N 2 e N then the operator AoT is of order | mi ~ 1,Tn ?~ 1 TO3 ~ ^ 2 | and class 
'cL iy 1 ) ' ^°^ e ^ na ^ the adjoint operator (A o T)* exists if iVi > max (mi, d 2 , di). 

2.3 System equations 

In the following we denote more shortly by A = (^j^ A^) a ^yP^ ca ^ e ^ emen ^ = ^ q 

of V. _ 

Let V\, vn and W\, w m be indeterminates in C°°(G x A) and in C°°(dG x A), respectively. 
Furthermore, let (Auj), (A 2 ik), (A3ij), {Auk) be operator matrices whereby are of type r + j4y + 
-By, A 2 ik are of type -Kifc, A^ij are of type Ty and Auk are of type Qik- We consider the control 
system 

N m 

^AujVj + ^AukWu = 0, i=l,...,iVi, (20) 
j=i k=i 

N m 

y^^j^ + A m w k = o, z = i, ..., 777,!. 

3=1 fe=l 

The state, input and output variables are not separated. 

The system ([20)) can be given in a matrix form as follows. Let (Auj) € M(N\ x N), (A 2 ik) £ 
M(N\ x m), (^4-3jj ) £ M(mi x iV) and (Aizfe) £ M(mi x m) be the matrices corresponding to the 



system (I2U1) . Denote tj = 



A 









(wi\ 


where d = 




... 













and denote 



(It;'! ;t:0 eM(( " i+m,)x(,v+m,)- <21) 



Then the system (I2U1) is equivalent to the equation 

Am = 0. (22) 

The operator A is a linear operator TZ N < m — > TZ Nl ' mi . We say the the operator A is in T>{N\ + 
mi, iV + 777) when it is of the form (|21|) . This notation uniquelly indicates the types of submatrices. 
Using the standard Frechet space topologies in C°°(G x A) and in C°°(dG x A) one has 
Assertion 4. The linear operator A : H N > m — > j£ N i> m i j s continuous. 

Example 3. In i/iis example we outline how to choose the above operators in the case of second 
order linear PDEs. Consider a linear coupled system of partial differential equations 

N 

A ij '(as, t, d)v 3 ■= 0, (23) 

3=1 
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for i = 1, ...,Ni, where Aij(x, t, d) are, for example, the second order operators 



A ij (x,t,d)= 

\a\<2 



where 



d 



d d d 

dxi ' '"' dx n ' dt 



(d x ,d t ). 



Denote Auj {x,t,d) = A. t j (x,t,d). 

We assume that the solution satisfies the following homogeneous boundary conditions 



N 



y^dij(x,t,d x )x 

3 = 1 



\dGxA = 0, 1 = 1, ...,mi 



(24) 



where di 3 {x, t, d x ) are first order partial differential operators 



i) 



dij(x,t,d x ) = ^ d i 3 k{x,t)— +d ljQ (x,t). 



k=i 



In this example we choose the operator matrix (Auj) such that 

( r+A u (x,t,d) ••• r+A 1N (x,t,d))\ 



(•Aitf) = 



\r+A Nll (x,t,d)) ••• r+A NlN (x,t,d))J 



where r + f := /|gxA is (as above) the restriction operator on G x A. The boundary operators 
Bij, K ikl Qik are zero operators and the operators A$ij corresponding to the boundary operators 
Tij are defined by 

/ r'd n (x,t, d x ) ■■■ r'd 1N (x,t, d x ) \ 



\r'd mi i(x,t,d x ) ■■■ r'd miN (x,t,d x )J 
where r' f := f\dGxA is the restriction operator on dG x A. As a conclusion we find that 



(25) 



A = 



{A Uj ) 
{M 3 ) 



(26) 



Remark 1. When we choose dij(x,t,d x ) = 0, I = 1, ...,mi, j = 1,...,N the function Vj does not 
satisfy any boundary condition on dG. 



Example 4. Consider the delay system 



dv 



— = A v + A k v{t - h k ) + B V 



(27) 



fc=i 



where < h x < ■ ■ ■ < h p , v(t) G C N , V(t) G C« and A Q , A k G M(N x JV), B a G M(N x q). 



Noting that 



f(t - h k ) = (27T)" 1 / fi^-^drj 

JR 

= (27T)- 1 / e-^fMe^dr, = P k f(t) 
Jr. 



(28) 
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we see that the system |<?7[) can be written in the form 



A ( y j = (29) 

where 

A=(|-A -ELi^ - fl o). (30) 
-Here is a pseudodiffererential operator with symbol Pk(v) = e~ lhk ' n . 

3 Parametrization of a control system applying algebraic 
methods 

3.1 Parametrizability of a control system 

Here we consider parametrizability of the above defined boundary value control systems. The 
control system 

N m 

AujVj + AgkWk = 0, i= 1, ...,Ni, (31) 

3=1 fc=l 
JV m 

^^3ij«j +5Z^ 4ZfcU7fe = 5 i = 1, mi, 
i=i i=i 

or more simply 

Au = (32) 
is said to be parametrizable, if there exists iV', m' £ N and a linear operator 



S = (if^l ^"T) e V{N + m,N' + ml) (33) 

\ \p-ikp ) {04kq)J 

such that 



Am = O u = S/ (34) 

(the components / m of / are X>-linear independent). The condition (|34[) means that the equation 
Am = determinates all compatibility conditions of the equation Sf = u. 

Example 5. Let G C R 2 . Denote d\ = j%£> d 2 — gfj and let I, I' be the identity mappings on 
C°°(G) and C°°{dG), respectively. Consider a system 

d 2 v\ - d\Vi - d 2 v 2 + v 2 = 0. 

The system can be put into the form (the inclusion of w\ is due to the notational convenience) 



d 2 -d x -d 2 + I 




Vl 

V 2 | : 0. 



The system has a parametrization ([201) 

'vA /d 2 -i o 

v i I = c?2 — di 

, Wl \ I 1 
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Example 6. Consider the control system related to Example\]\ 



L(D)vi = q(x)v 2 (36) 
d xVl (0) + di«i(0) = 0, a x «i(l) + da»i(l) = 

V = vi\dG- (37) 



This system can be put into the form 



[ T V I [ v 2 ) = 0. (38) 










-( 







Here I (resp. I') is the identity mapping on C°°(G) (resp. C°°(dG)). T is the trace operator 

\d x vi(x) + d2Vi(x), x = l 
Due to Example [JJ the system i38\) has a parametrization given by 

(P + P ) (<?(•)) \ 

I v 2 (40) 

r>(P + R )( q (-))J 

The homological algebra gives tools to study structurally the parametrizability. We describe 
the basic idea as follows. Suppose that T> is an algebra of suitable operators. Let Pi, P 2 , M be 
P-modules and let d\ : P 2 — > P\ and ■ Pi — * M be I?-homomorphisms. Recall that the sequence 
of modules 

P 2 — -*■ Pi M ^ 

is complex if im d\ C ker d 2 and it is exact if ker d 2 = im d\. The exactness means that d 2 v = 
if and only if v = dif and so the solutions of system d 2 v = can be parametrized by d\. The 
homology of the above complex is the quotient ker c^/im d±. The homology measures how much 
the sequence differs from being exact. Furthermore, using certain homology groups, Ext n (M, A), 
one can eventually specify how far a given module is e.g. from being projective. Using the groups 
Tor ra (M, A) one is able to measure how far a given module is from being flat. 



3.2 System modules 

As we found above, T> is an algebra of operators. Hence it is a (left) £>-module (for basic concepts 
of homological algebra see e.g. J22J, [TUl [13 E] ) . The module V has a unit 

x-(7 °) m 

where I is the identity operator C°°(G x A) — > C°°(G x A) and I' is the identity operator 
C°°(<9G xA)-> C°°(<9G x A). Hence V is a unitary P-module. 

We at first choose a subalgebra T>' of T> which contains an unit /. Furthermore, let T> lN be the 
direct product of modules 

V' N = V x • • • x V. (42) 
Then V' N is a (left) free X>'-module whose canonical basis is 

E x = {1,0,- ■■ ,0), ... E N = ((),••■ ,0,1). (43) 

Example 7. TTie following sets of operators are subalgebras of V 
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1. V = { [ r A ( D ) Q \ I A{D) is a PDO with constant coefficients}, 



2. V = { 



r+A{x,t,D) 




A(x,t,D) is a PDO with real analytic coefficients}, 



3. V = { ^dix't'l)) al' ' ' U ' r ' ° U ' h '''' D) PD()s C " (G " A '~ r,n ^"'""'^ " " 
R}, 



4- T>' = T>\, where T>\ is as in subsection \2.2.2i 
Suppose that the control system 

N rn 

AujVj + ^ A2ikW k 



0, i = l,...,Nt, 



(44) 



3=1 

N 



fc=l 



/] AsijVj + y j AukWk = 0, I = 1, ...,mi, 



fc=i 



is given. Let iV = maxjiV, m}, iVi = max{./Vi, mi}. The system (|44[) can be expressed equivalently 
in the form 

Cu ■ ■ ■ £, 



-in 



/•_ . . . r , 

• N 1 1 *-NiNj 



Wl 



= 0. 



(45) 



For example, suppose that N > m, Ni > mi. Define 

for 1 < j < to, 1 < i < mi, 



Cij — < 



■Auj 




Aiij 






:) 


















°) 





0/ 



for to + 1 < j < JV", 1 < i < mi, 
for 1 < j ' < m, mi + 1 < i < N\, 
for m + 1 < j < N, toi + 1 < % < N%. 



(46) 



Similarly we can reformulate the other cases. 
Example 8. The system 



/Am 


A 112 


Axxs 


A 2 11 


A 


Au = A121 


A\22 


-4l23 


Al2\ 


A 


\A 3n 


Az\2 


-4313 


A411 


A 



212 

222 
412 , 



V2 
V3 
Wl 

\w 2 J 



can be put into the form 



Wl 
V3 



= 0, 
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where 



C = 



£11 £12 £13 
£21 £22 £23 



( 




A 2 u\ 


Mm 


^212 






Am) 


1^312 


A412 




M121 


A221 \ 


Miaa 


-4.222 




v 


J 


I 







Let T>'(m x n) be the set of m x n- matrices whose elements are in I?'. Denote £ — (£y) G 



P'(iVi x TV) and denote uj 



w , 



, j = l, ...,N. Let it 



Then the system (J45j) is 



iV/ 



Cu = 0. 



Define a mapping .£ : T)' Nl — ► D w by 
.£D = £>£ = (Di 



Al 



"1JV 



(47) 



(48) 



"Nil ^NiNJ 

for D = (£>i • ■ ■ ^jVi) ^ T>' Nl . Then we find that X is a (left) P'-homomorphism that is, 

.L{L X D + L 2 D') = L^CD + L 2 .CD' 

for Li, L 2 S £>', D, D' G P'^ 1 . The image im(X) is exactly P /1? i£ and it is a submodule of V m . 
Hence we are able to define the (left) factor module 

M = coker(X) = V' w /V' Wl C. (49) 

The elements of M are of the form [D] = D + T>' Nl C Let ir : V N — > M be the canonical surjection 

■kD = [D]. (50) 

7r is a surjective P'-homomorphism and ker ir = T> lNl C and so we have an exact sequence 



V 



.IN 



M ■ 



0. 



(51) 



The module M corresponds to the control system (|4"4"]l in a sense that 



N 



J2^ij{Ej} = 0, l<i<N x , (52) 
3=1 

that is, the system equations are always valid in M for [Ej]. 

Remark 2. In £/ie case where the algebra P' is commutative we are also able to define a module 
M' := V' Wl /CV' 77 . In this case we find that 

C[D] = 0, VP e M' (53) 

that is, the system equations are valid in M' . Furthermore, we have an exact sequence 



V' 



IN 



D'Ni 



M' 







(54) 



A simple example of the commutative subalgebra is the first case in Example^ The commutativity 
is very restrictive property for P' . The module M' can be also defined in the case where T>' is only 
an Ore algebra. In this case |53) is valid. We do not consider module M' here. 
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The module M is finitely generated because for every element 
[D] = [{D 1 ••• %)]eM 

_ N 

[D] = (D 1 ••■ D w )+V'^C = J2D j [E j }. (55) 

The elements [Pj] do not, however, form necessarily a basis that is, every element [D] G M cannot 

be uniquelly expressed in the form [D] = Y^,jLi Dj[Ej]. By (fSTj) the module M is finitely presented. 

A P'-module M is free if it has a basis [Fx], [P r ]. Free module is always P'-isomorphic to 
the module P /r . A module M is projective free if there exists a P'-module P such that the direct 
sum M © P = F is free. The exact sequence of projective P'-modules Pj 

... ^p n _^vp n _ x ^ ... .P^PjAm -0 (56) 

is a projective resolution of M. The projective resolution always exists that is, we can choose 
projective P'-modules Po, Pi, P2, ... and P'-homomorphisms do, d\, g?2, •■■ such that 

do : Po — > M is onto 
dx : Pi — > ker do is onto 
c?2 : P2 — > ker c?i is onto 



and so on. 



3.3 Algebraic criterions for parametrizability 

Let S : R N ,m — > R N ' m be an operator as in section ETT1 



!?H eP(iV + m,iV' + m') 



(57) 



Similarly to the construction of C dehne the matrix S corresponding to S 



/Sxi ■■■ S l7 _ _ 

5= |e V'{N x AT'). (58) 

\S Wl ■ ■ ■ S m 

The homomorphism .S : V' N — > P /Af is also defined as X. 

Let A be a left P'- module. In this algebraic setting we say, more generally, that the system 
CX = 0, X G ^4^, is V -parametrizable in A N if there exists a matrix S given by l|58p such that 

CX = <s> X = (Xl, = SX' for X' G A 17 '. (59) 

In our case A may be 1Z or P', for example. When A = 7?. we get the parametrization concept in 
the sense of subsection 13.11 

Denote Hom(P, A) = {d\ d : P — > A is P' — homomorphism} for P'-modules P and A 
Hom(P, A) is a group with respect to addition but since V is not necessarily commutative Hom(P, A) 
is not generally a P'-module. We find that the mapping Hom(.£, A) : Hom(P /Ar , A) — > Hom(P /Ari , A) 
defined by 

Hom(.£, A)0 = o X (60) 

is a P'-homomorphism. 
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Lemma 1. The homomorphism Hom(X, A) can be calculated by 

JVi / JV \ 

[Hom(X,A)0](U) =X) A 53 ^3<KA) ( 61 ) 
»=i \3=i / 

forD = (D u ...,D 7fi )€V'^. 

Proof. For any D = Ysj=i DjEj G £> /Ar and for any </> S Hom(T>' N , A) one has the expression 

JV 

3=1 

and so 

[Hom(X, A)<j)](D) = {<j) o X)(D) = <t>(.CD) = <t>(DC) (62) 

N JVi / IV \ 

- 5] (D£), ^ao = 53 a 53 £ ^ ^ ) ( 63 ) 

3=1 i=l \3=1 / 

for L> e P'^ 1 as desired. □ 

Theorem 1. The parametrizability condition 

CX = 0<=>X = SX' for X' e A 77 '. (64) 

is equivalent to the equality 

ker Hom(X, A) = im Hom(.<S, A). (65) 
Proof. A. Suppose that |64|) holds. Then by Lemma Q] and by (|64|) we find that 

JV 

cj) e ker Hom(X, A) & Hom(X, A)<j> = ^ 53 L ^ E i) = 

3=1 

/0(A) \ F' 
«H i =5X'^^) = 5^.S jfe X^ 

wav)/ 

jv jv n' 

# 4>{D) = 5] DiftEj) = 53 A 53 (66) 

3 = 1 3=1 fc=l 

Define ip x > G Hom(P /17 ',^) by Vx'(-D) := EfeCi A X fc- Then we have X fe = V>X'(A) and so by 
Lemma [1] 

Iv jv' 

<t>{D) = 53 D^Sjrtx'iEk) = [Rom(.S,A)^ x ,}(D) & cj> = Rom(.S,A)^ x ,. (67) 

i=i k=i 



Since for any ip S Hom(P , A) 



'/'(A) 

r r v . (or A"' = j j | ((iS) 

vV-(Av') 
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the first part of the assertion is proved. 

B. Conversely, suppose that (|65p holds. Define <px(D) — 53j=i DjXj. Then by Lemma 1 

jV N 

CX = <s> C v X J = °' Vi £ y<M-Ej) = 0, Vi (69) 

3=1 3=1 

I, s 

^^A^Ai'M^j) = 0, VL> £ P'^ 1 <^Hom(.£,A)(fo: = 0. 

i=i j=i 

By (|65|) the equivalence (|69p is valid if and only if 

N W 

X = Hom(.S, & <f>x (D) = [Hom(.S, A)i/>] (£>) = ^ A\E 5 i*VW> V£> e ^ 

3=1 fc=l 

IV' / ^l) \ 

X 3 ; = ]T S jk ip(E k ) ^ X = SX', for X' = : . (70) 

Due to (f6"5]l this completes the proof. □ 
The homomorphism .S : V' N — » V' N is similarly defined as .C. Let N be the ©'-module 

N = V'~ W '/V ,W S. (71) 

For the definition and for the basic properties of the homology groups Ext ra (M, A) and Tor"(M, A) 

(see e.g. [HEBE]) 

We have 

Theorem 2. Suppose that there exists a matrix H58\) such that 

DC = D' & D'S = 0. (72) 
Then the parametrizability condition H65\) is valid if and only if 

Ext 1 (TV, A) = 0. (73) 
Proof. One sees that the assumption (j?2")) is equivalent to 

.CD = D' .SD' = 0. (74) 
Hence assumption (l?2"j) is equivalent to the exactness of the following sequence 

>- p/JV! — ^ V 'N v ,n' >. ft ^ o. (75) 

The resolution l|75|) implies the truncated dual sequence 

Hom(.S,A) _ Hom(.£,A) — , N 

>-Rom{V' N ,A) ^Eom(V' N ,A) >■ Uom(V Nl ,A) »- • • ■ (76) 

Ext 1 (N, A) :— ker Hom(.£, A)/im Hom(.5, A) — if and only if the condition ([65]) is valid. Hence 
Theorem Q] implies the assertion. □ 
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The condition (|73j) can also be studied without any explicit calculations. For example, since the 
module N is finitely presented, the condition f73|) is valid when the module N is flat. In the case 
where A is injective Ext 1 (iV, A) — 0. Especially, Ext 1 (TV, A) = in the case where N is projective. 
Note that neither V nor 1Z are generally injective modules. 

The conditions like (|73|) are interesting because they are independent of the used (projective) 
free resolutions. So they give valuable information for the instrinsic structure of the system. 

Example 9. Let G =]0, 1[ and A is any open interval o/R . Consider the system 

d 2 v\ d 2 V\ 
dt 2 dx 2 
«i(0,t) = ui(l,t)=0. 

In practice V\ may be the state variable and v-i may be the control variable. We have N = 2, N\ = 
1, mi — 1 and dG = {0, 1}. For notational convenience we may assume that m = 1 although no 
boundary value functions ui\ exist. Hence N = 2, N\ = 1. Denote d\ = d 2 = Jj. 
We see that 

c:= { c, &) -((^«-* s) s)). 

VFe seefc a possible operator S for which the assumption {7°2fy holds. Let D = I ^ ^ J G I? 



and (2T, 25") € P 2 ^ere D' = (^j , 23" = A }, 



A' 3 A'J 1 A'l 



We find that 



DC = D' (78) 



if and only if 



A 1 r+(d 2 -d 2 ) + A 2 r'I Q\ (-Ax 0\\_{{A[ A' 2 \ fA'{ A[ 

3 -^4/ \"^3 "^4 



A 3 r + {d 2 -d 2 ) + Air'I \-A 3 -\\A' 3 A'J U3 4" 



(79) 



i/ and on/y i/ 



^2 = >*4 = -4 = ^4 = 0, -Ai = A'{, -A 3 = A 3 

-A'{r + {dt - d 2 ) + A 2 r'I = A[, -Alr + (dl - d 2 ) + A 4 r'I = A 3 . (80) 
Let K be the potential operator 

Kg{x, t) = xg(l,t) + (1 - x)g(0, t). (81) 
Then we observe that r' Kg — g or r'K = I'. Multiplying the equation 

- A'{r+ (d 2 - d 2 ) + A 2 r'I = A[ (82) 

by K we get 

A 2 = A[K + A / lr+{di-d 2 )K. (83) 
Similarly we get from the equation 

- {d 2 - d 2 ) + A^' I = A! 3 (84) 

that 

Ai = A* 3 K + A^r+ (dl - d 2 )K. (85) 

Hence by {EU\) 

- A'{r + {d 2 - 8 2 ) + (A[K + A'{r+ (8 2 - d 2 )K)r'I = A[ (86) 
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and 

^r+{dl - df) + {A' 3 K + A^r + (d 2 2 di)K)r'l = A',. 
Expressing the equations [8b]) . |&7[ ) and 



(87) 



in the matrix form we have 



where 



A[ A' 2 



A'{ A 
A'i A'l 



h.o - Kr'I 
I 
(r+(di-df){-I + Kr'I) 




= 
















(0 





(89) 



(90) 



In section\K3\ we shall find that S is the V -parametrization in any module A 2 . 
3.4 Application of adjoints 

The adjoints can be applied also in the algebraic analysis which we describe in the sequel. In this 
subsection we assume that the required formal adjoints exist in the subalgebra T>' . 

Remark 3. If necessary we can apply order reducing to quarantee the existence of formal adjoints. 
For example, in the case where A = 1Z this can be described as follows. 

Suppose that T and T 1 are isomorphisms 1Z N — > 1Z N and lZ Nl —> lZ Nl , respectively. Then the 
control system Cu = is equivalent to the system 



(T o C o T)v! = 



(91) 



where Tu' — u. Especially, the operators T and T can be chosen to be order reducing that is, 
the orders and/or classes of the composed operators are less than the original operators (section 
\2.2.3\) . In parametrization we are able to use order reducing operators to quarantee the existence 
of adjoints. This is based on the following observation 



£!i = 0«-(r'o£o T)u' = u = S'f ^u=(To S')f &u = Sf 
where S' is the parametrization of the system (7~' o C o T)u' = and S = T o S 1 . 
Let L* be the formal adjoint of C 



(92) 



JVil 



IVilvV 



E V'(N 1 x N) 



(93) 



IN 



defined as in 12.2.21 The homomorphism .£ 
©'-module 

N = V' Nl /-ryA 



/V' N c*. 

Suppose that N has a free resolution with finitely generated modules 



V' Nl is similarly defined as Let N be the 

(94) 



V 



,/N 



•<S* ITT ■£* ITT 

*- jyiN q~)iNi 



N 



0. 



(95) 



A P'-module E is stably free if one can find a finitely generated free module F such that E © F 
is a finitely generated free module that is, E © F is isomorphic to 2?'™ for some n. A projective 
stably free module N always admits the free resolution ([95]) with finitely generated modules, for 
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example (pT]). In addition, in the case where V is Noetherian N has the free resolution with 
finitely generated modules (|95|) . For example, the subalgebra 

V = {{^ +A ^ ^ | A(D) is a PDO with constant coefficients}. 

is a commutative Noetherian integral domain. Besides this kind of simple cases, substantial further 
work must be done to analyze algebraic properties of more general systems. 

Remark 4. Suppose that A = V . Since CD = <^> D = SD' is equivalent to 

D*C* = D* = D'*S* (96) 

the parametrizability condition i59\) in the case A = T>' is equivalent to the existence of free reso- 
lution H95\). 



Suppose that the free resolution (|95|) exists. Let S := S** . As above we see that the condition 
Ext 1 (AT, V) = is equivalent to compatibility condition 

C*D = & & S*D' = 0. (97) 

Let D = D*, D' = D'*. Using the relation (£*£))* = DC = .CD and similarly (S*D')* = .SD' 
one sees that the compatibility condition (|9"7| is equivalent to 

.CD = D' <s> .SD' = 0. (98) 

Hence the condition 

Ext 1 (N, V) = (99) 
is equivalent to the assumption (|72|) of Theorem [2l 
Example 10. A. Let n = 2. Denote d\ = d 2 — gf^-. We choose 

V' = I ^ qJ I P{di,d 2 ) is PDO with constant coefficients^ . (100) 

Consider the following system (without boundary values) 

dxv 2 - d 2 v 1 = 0. (101) 

In this case we have 

c = (c u (*).«. -(£).«.-(£) 

where C\\ = jj J , £12 = ( q 2 jj J . W^e /mc? i/iai N = 2, N± = 1. 



5. T/ie formal adjoint exists and it is 



C = (^] CUK!) 



where C*i = ( ^ ^ J , £* 2 = f ^ 2 q j ■ TTie module N is given by 



N = V'/V' Z C*. (104) 
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Consider the existence of the resolution h9b^) . Let (D',D") £ T>' 2 . Denote D' 
where P' = P'{d u d 2 ), P" = P"(d 1 ,d 2 ). We find that 

X*(D',D") = 0<s> [D> D")C*=i) 



P' 0^ 
0, 



P" 




D" = 



(105) 



if and only if 



P' 0\ f-di 0\ fP" 0\ (b% 



0/ V 



o oy V o o 



= -P'di + P"d 2 = 0. 



Furthermore, the symbol of —P'd\ + P" d 2 is the polynomial —P' (£i£2)£i + P" (^1^2)^2 a71 <i one 
knows from algebra that there exists a polynomial Q such that (Ulf . P- 113) 

-P'(a6)a + J 3 "(a,6)6 = o^P'(66) = Q(ei,6)6, i y/ (6,6) = Q(fc,6)&. (106) 

Hence the condition —P'd\ + P"d 2 — is equivalent that 

P' = Q(di,d 2 )d 2 , P" =Q(d 1 ,d 2 )d 1 . 



(107) 



Denote S* = 



d 2 0\ /9 X 











^4s a conclusion we see that 



.£*(£>', £>") = ^ (D , D ) = Q{dt,d 2 ) 
(D\ D") = .S*Q. 



d 2 0\ /9i 











Q(dud 2 )S* 



(108) 



Hence N = 1 and £/ie /ree resolution H95\) with finitely generated modules exists. 
C. It is easy to see that Ext^iV, T>') = 0. Furthermore, 



Ext 1 (TV, ft) = if and only if 



where U = 
and only if 



Vi 
Wi 



£ 1Z and u\ 



s = 


S** = - 


su 


"GO 












, u 2 = 











°)\ 













o\ 


V 


\ 


V) 



VVi = (wi, U2) 9iu 2 - 9 2 «i = 
which is valid if the set G is convex, for example. 



(109) 



Ul )=0 (110) 
u 2 J 

£ 1Z. Hence we see that Ext 1 (TV, TZ) = if 

(111) 



3.5 Projectivity 

We finally shortly treat the projectiveness of the underlying modules. The module M is projective 
if and only if Ext 1 (M, A) — for any left module A. The projective dimension is defined by 

P - dim M = inf {Ext" +1 (Af, A) = for any A}. (112) 

n>0 

We see that P-dim M = if and only if M is projective. For non-projective modules M projective 
dimension measures how far M is from being projective. 

The following theorem gives a tool to show the projectiveness in our case. 
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Theorem 3. The module M is projective if and only if there exists a T>' -homomorphism .V : 
v fN v ,7i l ( g0 calM m j guch that 

X o .V o X = X. (113) 

The proof of theorem can be found e.g. in 20J . The condition (|113j) is equivalent to the 
existence of a matrix V G T>'(N x N\) such that 

CVC = C. (114) 

Corollary 1. Suppose that there exists a T>' -homomorphism .V : T>' Nl — > D' N (so called left 
inverse ) such that 

.PoX=l Wi (115) 
where Xjy is the identity in T>' Nl . Then the module M is projective. 

The condition (| 1 1 5(1 is equivalent to the existence of a matrix V G V(N x N\) such that 

It o---o \ 

CV = X 7 f 1 :=\\ : : : e V'(N 1 x N x ). (116) 
\0 • • • ij 

Projectivity gives the following sufficient condition for the parametrizability 
Theorem 4. Suppose that there exists a matrix V G T)'(N x N\) such that 

CVC = C. (117) 

Then the control system CX = is T>' -parametrizable in A N . 
Proof. From (|117[) we get 

- VC) = 0. (118) 

Hence the operator S = — 7 ? £ satisfies the following condition: If X = SX' then CX = 0. 
Conversely, suppose that CX = 0. Then 

X = (% - 7 7 £)X + (VC)X = SX (119) 

which completes the proof. □ 

A special case of the Theorem |4] is 

Corollary 2. Suppose that the matrix C has a right inverse that is, there exists a matrix V G 
V(N x Ni) such that 

CV=l 77i . (120) 
Then the control system CX — is T>' -parametrizable in A N . 

Under the assumptions of Theorem [4] the parametrization is given by 

SX' = (% - VC)X\ X' G A W . (121) 

We find that the projectiveness is a structural property under which the parametrizability does 
not depend on the module A. 

Remark 5. A. The module M is flat if and only if Tor-functor is exact. Projective module is flat 
but there exist nonprojective flat modules. Flatness is a useful additional concept in the study of 
these system modules. 

B. For example, for integral domains T>' one sees by a direct computation that &98\) or |7£| ) is 
equivalent to the torsion freeness of M that is, D[D] = 0, D G V , [D] G M if and only if D = 
or [D] = 0. For nonintegral domains this need not be valid. Recall that 

free C stably free C projective C flat G torsion free. (122) 
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Example 11. Reconsider the Example^ Recall that 

The corresponding complex related to the module M is 

*■ M 



■V—^-V 2 



-I 




0. 



(123) 



(124) 



Furthermore, D.C = for D 



Ai A 



£ T> if and only if 



Ai = A 3 = A 2 r'I = A 4 r'I = 0. (125) 
The last two equations imply that also A% = A± = and so D = 0. Hence we get the free resolution 

.C 



0- 



The truncated dual sequence is 



0- 



Hom(D 2 ,2?) 



M ■ 



Hom(X>, V) 



•0. 







(126) 
(127) 



where $((/>) = Hom(.£,X>)0 = 0o X. TTie homomorphism $ is surjective: Due to Lemma\]\we see 
that <&(</>) = "0 /or "0 6 Hom(2?, £>) j/ and onfo/ if 

2 

5^£^(B i )=^(Bi). (128) 

3=1 

Choose (j){D 1 ,D 2 ) = (D u D 2 )Vi/j(Ei) where 



V 



( (0 K 

\Q 

f-I r + (d%-df)K 

Wo o 

Here, as above, K is the potential operator 



\ 



S = 



(130) 



Kg(x, t) = xg(l, t) + (l- x)g(0, t) (129) 

(observe that r'Kg = g). Then il28\) holds and so $ is surjective. Hence Ext 1 (M, T>) = 0. 

Actually M is projective. The sequence 1126( 1 splits and the lift can be chosen to be .V . We can 
compute the parametrization by H121)) and the result is 

I') lo 0) 
(0 0\ (I o\ 1 L 
\{o o) [o I' J J 

( (I - Kr'I 

^ o r 

f(r+(di-df){-I + Kr'I) 

VV o o 

By U30\) the parametrization in the sense of the subsection lS. 1\ is given elementwise (that is in the 
case A = 7i) by 

vi = }\-Kr'h 















G 





V2 



{d 2 2 - ODh - r+{d 2 2 - d\)Kr'h, fx e C°°(G x A) 
g u g\ e C°°(dG x A). 



(131) 

Note that Kr' f\ = xfi(l, •) + (1 — x)fi(0, •). The last equation of M31]) is due to the notational 
covensions and it is superflous. 
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Example 12. Let n = 2 and let Av% = + Consider the system 



(1 - A)vi = (132) 
dvi 

— | aG = W i. 

(133) 



Here we have N — m = N± = mi = 1 and so N = N\ = 1. In addition, 



where I' is the identity operator on C°°(<9G). Note that = r '( l/1 a§r ^ 2 3§j)' ^ e see ^ a 
matrix PeP such that 

CV=Tt = (i °,Y (135) 
Let V = _ -T-) • /md £/ia£ i/ie condition \135}) is equivalent to 

d d 

(1 - A)Vi =1, (1 - A)V 2 = 0, r'—Vx -V 3 = 0, r'—P 2 -V± = I'. (136) 

ov ov 

We can choose V 2 — and then V4 = — i 7 . ie£ Pi = r + A where r + A is the pseudo- differential 
operator with symbol 1/(1 + ||£|| 2 ). Then (l-A)Vx =1. Fromr' -^P 1 -P 3 = 0wegetT 3 = r'-^'P 1 . 
Hence we can choose 

*- -"')• (137> 

XTie corresponding parametrization is 

*-*-«.(ji)-(,^i)e^i) <«> 



7-(r + A)(l-A) 
-r'|(r+4)(l-A) + r'|; 

Elementwise the parametrization is given by 



(139) 



wi = (/-(r + A)(l-A))/i (140) 
«, 1 = _r'|-(r + A)(l-A)/i + r / ^. 

The matrix V can be chosen several other ways. For example, we can choose V\ such that 
(1 — A)Vi = I, r' -§ifPi = 0. Then V\ = r + A + B where B is a singular Green operator (JBj). 
Choosing V 2 = we get 

V=('*\ +B ",)■ (141) 



-I' 

Parametrization can be computed and it is of the form 



wi = (/-(r+A + B)(l-A))/i (142) 

wi=r — — . 
Ov 
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